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Summary. The fate of the cosmological singularity in the Taub model is discussed within 
the two frameworks. An internal time variable is ruled out and the only remaining degree 
of freedom (the anisotropy) of the Universe is quantized according to such schemes. The 
resulting GUP Taub Universe is singularity-free, differently from the second case, where the 
classical singularity is not tamed by the polymer-loop quantum effects. 



1 Introduction 

Two different quantum cosmology approaches are applied to the Taub model. The study 
is performed at the classical and at the quantum level in both schemes. In particular, the 
generalized uncertainty principle (GUP) and the polymer (loop) frameworks are implemented 
to this system. In the first case [1] , the cosmological singularity appears to be probabilistically 
suppressed, while, in the second one [2], the Universe is still singular. Such a feature then 
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allows us to better understand the avoidance of the cosmological singularity in other different 
quantum gravity toy models. 

The Taub Universe arises as a particular case of the Bianchi IX model, i.e. the most 
general scheme allowed by the homogeneity constraint (for reviews sec [3]). It is obtained by 
restricting the dynamics to that of a one-dimensional particle bouncing against a wall, when 
only one degree of freedom (the Universe anisotropy) is taken into account. The relevance 
of the Taub Universe in quantum cosmology is then duo to the fact that it is a necessary 
step towards the Bianchi IX model, being a generalization of other isotropic models. In 
particular, it has been used to test the validity of the minisuperspace scheme [4] and to 
explore the application of the extrinsic cosmological time [5] . 

The paper is organized as follows. In Section 2 the Taub model and the formalisms are 
reviewed. Section 3 and 4 arc devoted to the classical and quantum analysis respectively. 
Comparisons with other approaches follow in Section 5. We adopt h = c = 16nG = 1 units. 



2 Taub Universe and formalisms 

The Taub cosmological model is a particular case of the Bianchi IX Universe, of which 
only one anisotropy is taken into account. The Bianchi IX model, together with Bianchi 
VIII, is the most general spatially-homogeneous model and is described by the line element 
ds^ = N^dt'^ - e^" (e^^)^ . to' ® u}^ [3]. Here N = N(t) is the lapse function, lj' = cuidx" arc 
the 5*0(3) left-invariant 1-forms, a = a{t) describes the isotropic expansion of the Universe 
and 7i-, = ^ij{t) is a traceless symmetric matrix, which determines the anisotropics via 
7±. The dynamics of this Universe towards the singularity is described by the motion of 
a two-dimensional particle (the two physical degree of freedom of the gravitational field) 
in a dynamically-closed domain [3]. In the Misncr picture, such a domain depends on the 
time variable a, while, in the Misner-Chitre one, it becomes stationary in time. Performing 
the ADM reduction of the dynamics (according to which the classical constraints are solved 
with respect to the given momenta before implementing any quantization algorithm), an 
effective Hamiltonian is obtained, which depends only on the physical degrees of freedom 
of the system. In particular, the scalar constraint is solved with respect to the momentum 
conjugated to the time variable pr (we adopt the time gauge f = 1) and, performing another 
change of variables, we obtain 

-Pr = Hi^M = vy/pl+Pv- (1) 

The dynamics of such a system is equivalent to a billiard ball on a Lobatchevsky plane and 
the three corners of the Misncr scheme are replaced by the points (0, 0), (—1, 0) and v ^ oo 
in the (tt, t;)-plane, as in Fig. 1. The Taub cosmological model is described by 7- = 0. The 
dynamics of this Universe is equivalent to the motion of a particle in a one-dimensional 
closed domain. Such a domain corresponds the choice of only one of the three equivalent 
potential walls of the Bianchi IX model. The ADM Hamiltonian (1) rewrites i^JiJM = ''JPv, 
where £ [1/2, oc). This Hamiltonian can be further simplified defining the new variable 
X = \nv, and becomes 

Hadm =Px=P (2) 

where the configuration variable x is related to the Universe anisotropy 7+ by the equatiou 
7+ = e^~^{e^^ — 3/4)/\/3. The Hamiltonian (2) will be the starting point of our analysis. It 
is worth noting that the classical singularity now appears for r — > 00. 
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Fig. 1 The dynamical- allowed domain in the {u, t;)-plane where the dynamics is restricted. 

2.1 GUP quantum mechanics 

Some issues and results of a non-relativistic quantum mechanics with non-zero minimal 
uncertainties in position are briefly reviewed [6] . In one dimension, we consider the Heisenberg 
algebra generated by q and p obeying the commutation relation 

[q,p] = i(l + (3) 

where /9 > is a deformation parameter. This commutation relation leads to the generalized 
uncertainty relation AqAp > [l + f3{Ap)'^ + P{p)'^) , which appears in string theory [7]. The 
canonical Heisenberg algebra can be recovered in the limit /3 = 0, and the generalization to 
more dimension is straightforward, leading naturally to a "noncommutative geometry" for 
the space coordinates. 

It is immediate to verify that a finite minimal uncertainty in position Aqmi„ = is 
predicted. The existence of a non-zero uncertainty in position is relevant since it implies 
that there cannot be any physical state that is a position eigenstate. In fact, an eigenstate 
of an observable necessarily has a vanishing uncertainty on it. Although it is possible to 
construct position eigenvectors, they axe only formal eigenvectors and not physical states. 
The deformed Heisenberg algebra (3) can be represented in the momentum space as 

pVb) = pV'b), = i(i + /3p^)5pV(p), (4) 

on a dense domain S of smooth functions. To recover information on position we have to 
study those states, wliieli realize the maximally-allowed localization. Such states IV'™')) which 
are proper physical states around a position ^, have the proprieties (V'™'klV'c"') — C S'Hd 
(Zig)|^mi) = Aq,nin- We can project an arbitrary state lip) on the maximally-localized states 

\ipc^') to obtain the probability amplitude for a particle being maximally localized around the 
position ^ (i.e. with standard deviation Aqmin)- We call these projections the quasiposition 
wave function tp{C) = (V'c"IV')i and explicitly we have the generalized Fouries transformation 

As (3 —^ 0, the ordinary position wave function ^{Q = (C|^) is recovered. 

2.2 Polymer quantum mechanics 

The polymer representation of quantum mechanics consists in defining abstract kets, labeled 
by a real number and assumed to form an orthonormal basis, and then considering a suitable 
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finite subset of them, whose Hilbert space is defined by the corresponding inner product [8]. 

This procedure helps one gain insight onto some particular features of quantum mechan- 
ics, when an underlying discrete structure is somehow hypothesized. The request that the 
Hamiltoniam associated to the system be of direct physical interpretation defines the polymer 
phase space, and the continuum limit can be recovered by the introduction of the concept 
of the scale [9]. 

In the particular case of a discrete position variable in the momentum polarization, the 
Hamiltonian variable p cannot be implemented as an operator, so that some restrictions on 
the model have to be required. If the set of kets is restricted by the introduction of a regular 
graph 'y^g , the kynetic term of the Hamiltonian is approximated by the polymer substitution 

p -> — sin(/iop), (6) 
Ho 

where the incremental ratio is evluated for an exponentiated operator. The Hamiltonian 
operator H^^^ , which lives in Ti.^/^^ , reads 

^^.o = £ + m (7) 

The definition of a scale, Cn, cables one to approximate continuous functions with functions 
that are constant on the intervals. As a result, at any given scale C„, the kinetic term of the 
Hamiltonian operator cam be approximated, and effective theories at given scales are related 
by coarse-graining maps. 



3 Deformed classical dynamics 

The ordinary Taub model can be interpreted as a massless scalar relativisitic particle moving 
in the Lorentzian minisuperspace (t, a;)-plane, whose the classical trajectory is its light-cone. 
More precisely, the incoming particle (r < 0) bounces on the wall {x = xq — ln(l/2)) and 
falls into the classical cosmological singularity (t oo). Investigating the modification of 
the dynamics within the two frameworks will show that the two behaviors can be interpreted 
as complementary. 



3.1 GUP framework 

The GUP-classical dynamics is contained in the modified symplectic geometry arising from 
the classical limit of (3), as soon as the parameter (3 is regarded as an independent constant 
with respect h. It is then possible to replace the quantum-mechanical commutator (3) via 
its Poisson brackets, i.e. —i[q,p] =4> {q,p} = (1 -|- The Poisson brackets for any two- 

dimensional phase space function are 



dq dp 

Applying this scheme to the Hamiltonian (2), we immediately obtain the equations of motion 
for the model [1], 

x(t) = (1 + i3A'^)T + cost, p{t) = cost = A, (9) 
where x G [a;o,oo). Therefore, at the classical level, the effects of the deformed Heisenberg 
algebra (3) on the Taub Universe are as follows. The angular coefficient is (1 + f3A^) > 1 for 
/9 7^ 0, and thus the angle between the two straight lines x(t), for r < and r > 0, becomes 
smaller as the values of /3 grows. The trajectories of the particle (Universe), before and after 
the bounce on the potential wall aX x = xo = ln(l/2), are closer to each other then in the 
canonical case {/3 = 0). 
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3.2 Polymer framework 

The polymer-classical dynamics relies on the substitution (6) in the Hamiltonian of the model 
(2). This way, the equations of motion rewrite [2] 

X = {x, H} = cos(ap), p = {p, H} = 0, (10) 

where a dot denotes differentiation with respect to the time variable t. The equations of 
motion are immediately solved as 

x{t) = cos{ap)T, p{t) = A, (11) 

where A is a constant. In the discrotizcd (polymer) case, i.e. for a 7^ 0, the one-parameter 
family of trajectories flattens. In fact, the angle between the incoming trajectory and the 
outgoing one is greater than 7r/2, since p € (— 7r/o, tt/o). As these trajectories diverge rather 
than converge, we expect the polymer quantum effects to be reduced with respect to the 
classical case, as we will verify below. 



4 Deformed quantum dynamics 

The quantum dynamics of the Taub Universe is here investigated according to the two 
different approaches. Particular attention is paid to the wave-packet evolution and the con- 
sequential fate of the classical cosmological singularity. In both frameworks, the variable r 
is regarded as a time coordinate and therefore (t,Pt) are treated in the canonical way. The 
deformed quantization (GUP or polymer) is then implemented only to the submanifold de- 
scribing the only degree of freedom of the Universe, i.e. the phase space spanned by {x,p). 
We then deal with a Schrodinger-like equation 

id.^{T,p)^HloM9{r,p), (12) 

where the operator H'adm accounts for the modifications due to the two frameworks. We 
have to square the eigenvalue problem in order to correctly impose the boundary condition: 
in agreement with the analysis developed in [10], we make the well-grounded hypothesis that 
the eigenfunctions form be independent of the presence of the square root, since its removal 
implies the square of the eigenvalues only. The wave packets, which are superposition of the 
eigenfunctions if'(r, x) = dkA{k)tpk{x)e~'''''^ , are then constructed for both models, taking 
A{k) as a Gaussian-hke weighting function. The differences between the two approaches are 
due to the features of the eigenfunctions ipk{x). Analyzing such an evolution, we show that 
the GUP Taub Universe appears to be probabilistically singularity-free, differently from the 
polymer case, where the singularity is not tamed by the cut-off-scale effects. 



4.1 GUP framevirork 

We now analyze the model in the GUP approach [1]. As explained before, we lost all infor- 
mations on the position itself, so that the boundary conditions have to be imposed on the 

quasiposition wave function (5), i.e. 4'iCo) = (where ("0 — (V'™'la;-()|-0™'), in agreement with 
the previous discussion). The solution of the eigenvalue problem is the Dirac J-distribution 
V'fc(p) = ^iP^ ~ fc^)i and therefore the quasiposition wave function (5) reads 



MO = 



W4W^ [exp {ij^t.n-\^pk)) -exp (.(^tan-(^fe) 



, (13) 
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where A is a constant and the boundary condition V'(Co) = has been imposed. The deforma- 
tion parameter /3, i.e. the presence of a non-zero minimal uncertainty for the configuration 
variable, is responsible for the GUP effects on the dynamics. The physical interpretation 
of 13 is then a non-zero minimal uncertainty in the anisotropy of the Universe. To better 
understand the modifications induced by the deformed Heisenberg algebra on the canonical 
Universe dynamics, we have to analyze different /3-regions. In fact, when f3 becomes more and 
more important, i.e. when we are at some scale that allows us to appreciate the GUP effects, 
the evolution of the wave packets is different from the canonical case. More precisely, these 
effects are present when the product k^y/P becomes remarkable, i.e. when koy/P ~ 0(1), and 
therefore when (3 is comparable to l/k^. In fact, the correct semiclassical behaviors of the 
model far away from the singularity is described by wave packets peaked at energies much 
smaller then l/\fP [1]. In particular, for fco = 1, we can distinguish between three different 
/3-regimes: 

• /3 ~ 0(10"^) regime. The wave packets begin to spread and a constructive and de- 
structive interference between the incoming and outgoing wave appears. The probability 
amplitude to find the Universe is still peaked around the classical trajectory. 

• /3 ~ 0(10^^) regime. It is no more possible to distinguish an incoming or outgoing wave 
packet and, at this level, the notion of a wave packet following a classical trajectory 
becomes meaningless. 

• /3 ~ 0(1) regime. A dominant probability peak "near" the potential wall appears. There 
are also other small peaks for growing values of but they are widely suppressed for big- 
ger /3. The motion of wave packets shows a stationary behavior, i.e. these are independent 
of r. See Fig. 2. 

Following this picture we are able to learn the GUP modifications to the WDW wave packets 
evolution. In fact, from small to big values of j3, we can see how the wave packets escape 
from the classical trajectories and approach a stationary state close to the potential wall. 




Fig. 2 Wave packets |*'(r, C,)\ in the GUP framework as /3fco = 1 (fco = 1 and a = 4). 

Such a behavior is, in some sense, expected from a classical point of view. In fact, at 
classical level the ingoing and the outgoing trajectories shrink each other. So a quantum 
probability interference is a fortiori predicted. On the other hand, the stationarity feature 
exhibited by the Universe in the (/3 ~ 0(l))-region is a purely quantum GUP effect. Such a 
behavior cannot be inferred from a deformed classical analysis. From this point of view, the 
classical singularity (r oo) is widely probabilistically suppressed, because the probability 
to find the Universe is peaked just around the potential wall. This way we claim that the 
GUP-Taub Universe is singularity-free. 
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4.2 Polymer framework 

We now analyze the model in the polymer approach [2]. For the quantum analysis of the 
model, we choose a discretized x space, and solve the corresponding eigenvalue problem 
in the p polarization. Considering the time evolution for the wave function ^ as given by 
'I'k{p,'r) = e~^'''^4>k{p) and the results of [10], we obtain the following eigenvalue problem 



(p^ - k'^)Tpk{p) = 



2 

— (1 - cos(ap)) - 
a-' 



Mp), (14) 



solved by 



fe^ = fc^(a) = (1 - cos(ap)) < kl^ax = ^ (15a) 

i^k.aip) = AS{p-pk,a)+B5{p+Pk,a) (15b) 

ipk.a{x) = A [exp{ipk,ax) - exp{ipk,a{2xo ~ x))] : (15c) 

(15b) is the momentum wave function, with A and B two arbitrary integration constant, and 
(15c) is the coordinate wave function, where an integration constant has been eliminated by 
imposing suitable boundary conditions. Moreover, we have defined the modified dispersion 
relation 

1 / . X 

Pk.a = - arccos II — I (16) 

from (15a). Furthermore, we stress that k^ is bounded from above, as illustrated in (15a), 
but it is its square root, considered for its positive determination, which accounts for the 
time evolution of the wave function. 

We now construct suitable wave packets 'l'{x, r) taking into account the previous dis- 
cussion (note that a maximum energy kmax is now predicted). Three relevant cases can be 
distinguished: 




Fig. 3 The spread polymer wave packet |^'(a;,r)| as koa = 1/2 (a — 50, ko = 0.01, 
a = 0.125). 



• koa ~ 0(1) and peaked weighting function. The resulting wave packet is well approxi- 
mated by a purely monochromatic wave. A small interference phenomenon between the 
wave and the wall is then predicted. 

• koa ~ 0(1) and spread weighting function. A strong interference phenomenon between 
the wave and the wall now appears. Nevertheless, this interference phenomenon is not 
able to probabilistically tame the singularity, as it takes place in the 'outer' region, in 
a way complementary to that of the GUP approach (see Fig. 3). The polymer- Taub 
Universe is then still a singular cosmological model. 
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• koa <C C(l) regime. This can be considered as the semi-classical limit of the model. In 
fact, differently from the other cases, the value of ko around which the wave packet is 
peaked is not arbitrary, but constrained by the characteristic scale a under investigation. 
The ordinary WDW behavior is therefore recasted. 



5 Comparison with other approaches 

The Taub cosmological model offers a suitable scenario, where different quantization tech- 
niques can be applied. In fact, it is possible to single out a time variable, so that the anisotropy 
describes the real degree of freedom of the Universe. It is therefore reasonable to investigate 
the fate of the cosmological singularity without modifying the time variable. The comparison 
with analysis of the cosmological singularity in other cosmological models outlines how the 
features of the Taub model allow one to pick the cut off effects out of those due to the choice 
of the Hamiltonian variables. 

In the cosmological isotropic sector of GR, i.e. the FRW models, the singularity is re- 
moved by loop quantum effects. The wave function of the universe exhibits a non-singular 
behavior at the classical singularity, and the big-bang is replaced by a big-bounce, when a 
free scalar field is taken as the relational time [11]. The scale factor of the universe is directly 
quantized by the use of the polymer (loop) techniques, so that the evolution itself of the 
wave packet of the universe is deeply modified by such an approach. The Hamiltonian con- 
straint does not allow for a constant solution of the variable conjugated to the scale factor, 
so that it is not possible to choose a scale, such that the polymer modifications are negligible 
throughout the whole evolution, so that the comparison with the ordinary representation is 
not always possible. Anyhow, we stress that, for the Taub model, the cosmological singularity 
is probabilistically suppressed, regardless to the fact whether the system can appreciated or 
not the cut off during the whole evolution. 

In [12], all the degrees of freedom of the Bianchi cosmological models in the ADM reduc- 
tion of the dynamics axe quantized by loop techniques. In particular, also the time variable, 
i.e. the Universe volume, is treated at the same level as the others. In most cases, the time 
variable is defined by a phase space variable, i.e. it is an internal one. As a result, also the 
Bianchi Universes are singularity-free [13]. In this respect, our analysis is based on consider- 
ing the time variable as an ordinary Heisenberg variable, while the cut off is imposed on the 
anisotropy only. 

The GUP dynamics of other cosmological models has been investigated in different ap- 
proaches. In particular, the big-bang singularity appears to be tamed by GUP effects showing 
a stationary behavior of the wave packets [14] . Such a prediction is in agreement with those 
achieved in a noncommutative quantum cosmology [15]. However, in order to predict a big 
bounce a la LQC, a Snyder-deformed quantum cosmology has to be addressed [16]. As the 
last point, it is interesting to notice that the GUP-Mixmaster Universe is still a chaotic 
model [17], as opposite to the LQC one [18], the difference being essentially based on the 
application of the deformed scheme to the time variable too. 

Acknowledgments. The "Angelo Delia Riccia" Fellowship and the "Sapienza GUN 2" 
Fellowship are gratefully acknowledged. 
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